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ABSTRACT 


The hypersonic laminar two-dimensional flow over a stepwise- 
accelerated semi-infinite flat plate at zero angle of attack is analysed. 
The boundary layer governing momentum and energy equations are trans- 
formed in terms of the dimensionless stream function and the dimension- 
less total enthalpy (S). These coupled equations have been solved 
numerically for Prandtl number unity, pressure gradient parameter 
B= 0.286 and 0.4,0<S <1, and0 << 0.3 withe= (1+ X/V,T), 
The obtained solutions are utilized to compute the time-dependent dis- 
placement thickness which is subsequently used to obtain the strong 
interaction induced pressure (according to the Lighthill's Piston Theory) 
for 5 Me < 10. For the same range of Mach numbers, the transient 
contributions to the boundary-layer velocity, the boundary-layer temper- 


ature and to the shear stress at the wall are also given. 
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k thermal conductivity 
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m an integer in numerical method 
M Mach number 
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Pr Prandtl number 

q heat transfer 


Q count of successive solutions of the momentum equation in the 


procedure of solution 


R gas constant 

R coefficient, defined by eq. (A.31) 

Re Reynolds number 
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transformed time variables 


temperature 


transformed time variables 
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E modified similarity independent variable associated with 


time and space, defined by eq. (3.3) 


n similarity independent variable associated with space 

v dynamic viscosity coefficient 

v kinematic viscosity coefficient 

le) density 

1: shear stress 
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CHAPTER I 
INTRODUCTION AND DISCUSSION OF THE LITERATURE 


1.1 Statement of the Unsteady Hypersonic Flow Problem 

There are many situations where a quantitative description of 
the transient development of a boundary layer would be desirable and 
valuable; a few of these are aircraft and missiles in unsteady flight, 
oscillating wings, unsteady nozzle flow, blades rotating in non-uniform 
air streams, etc. In evaluating the forces acting on the components 
of aircraft and missiles during their maneuverability and control, 
viscous interactions have to be analysed. The viscous aspects of shock 
wave interaction phenomenon for such a case require the study of time- 
dependent boundary layers. The sharp-leading-edge bodies shall be 
considered for purposes of our investigation. 

The problem to be discussed here, considers the interaction 
between the viscous and inviscid effects, when the high Mach number 
flow over a flat plate at zero angle of attack is stepwise-accelerated 
by a small amount in the direction of the initial flow. The induced 
pressure due to the viscous interactions will be obtained from the 
analysis of the temporal velocity and enthalpy boundary layers. Ex- 
pressions will also be obtained for the transient shear-stress and heat 


transfer at the wall. 
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1.2 Review of the Associated Literature 

The first treatment of the transient boundary layer was pre- 
sented by Blasius [1] in 1908, who considered the "impulsive" start of 
bluff bodies such as cylinders. More recently, the review articles by 
Stewartson [2] and Rott [3] have covered the field in greater detail. 
Some of this literature, which has direct bearing to our work, will be 
discussed here. 

Moore [4] considered the case of compressible laminar flow 
over a semi-infinite flat plate. In his analysis he found a group of 
parameters, whose magnitude determines the nature of flow unsteadiness 
for an arbitrary velocity U(t). Ostrach [5] extended Moore's work of 
zero wall heat transfer to the case of an isothermal surface. These 
are the asymptotic solutions for large times. In a recent paper Yang 
and Huang [6] supplemented the above with the solutions for small times. 
The boundary layer solutions of references [4], [5] and [6] have been 
utilized by Zien and Reshotko [7] to obtain the transient weak viscous 
interaction pressures. As a first approximation to the induced surface 
pressure, they have taken it equal to the acoustic pressure on a one- 
dimensional piston moving at a variable, but low, speed. 

Stewartson [2,8], while considering the impulsive motion of a 
flat plate, gave an interesting explanation of how the initial Rayleigh 
solution, having no x-dependence, could possibly change over to the 
x-dependent but time-independent Blasius solution. He obtained two 


solutions, one for t(= Ut/x) > 1 and the other for t < 1 with a singu- 
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larity at t = 1. He employed two different methods, namely, the Oseen's 
approximation and the momentum integral method to obtain these solutions. 
Stewartson argued that for t << 1, since a disturbance from the leading 
edge has not arrived yet, Rayleigh-type flow prevails; whereas, at 
t = | the leading edge effect is felt suddenly and the flow acquires 
its x-dependence through an essential singularity at this point. The 
existence of singularity was further supported by Smith [9] who gener- 
alized Stewartson's work for a wedge. 
Lam and Crocco [10] investigated the boundary layer induced 

by a shock travelling down a semi-infinite flat plate. At time t = 0 
the shock is at the leading edge of the plate, x = 0. For the special 
case of a very weak shock the problem becomes identical to Stewartson's 
impulsive plate problem. Using the Prandtl boundary layer equations 
under the Crocco transformation they give sufficient conditions for 
unique solution to the governing equations for the case py = constant 
and demonstrate the feasibility of enforcing two streamwise boundary 
conditions for such a problem. The boundary layer is divided in two 
parts: 

(1) a region where the solution of Mirels [11] is valid (t < 1); 

(2) the region, t >], 
where both the leading edge effects and the effect of the shock solu- 
tion downstream must be considered. The results are directly analogous 
to Stewartson's. In the region t > 1 the governing equations, although 


parabolic in nature, admit boundary conditions which are usually as- 
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sociated with elliptic equations. This type of behavior has been termed 
“singular parabolic" by Gevrey [12] and also Lam and Crocco ElO}s- re 
is in this second region that the "elliptic" boundary conditions must 
be enforced; i.e., boundary conditions must be given on t = 1, ~ for 
0 < n(= 7 < 1 and onn = 0, 1 for 1 <t< ©, Though the complete 
solution ts the problem was not given by Lam and Crocco, the detailed 
discussion of the properties of the solution was provided. 

Stewartson's singularity was also encountered by Cheng [13] 
who investigated the problem of a sharp-edged flat plate starting from 
rest and following a power law motion. Cheng expands the stream function 


in powers of quantity 


p= Sint) with u(t) = at” . 
At : 


He obtains a solution which is valid for large times (i.e. & << 1). 

For large € the solution is of the Rayleigh type. Later Cheng and 
Elliott [14] extended this to arbitrary plate velocity if the plate was 
started initially from rest. 

Using the approximate momentum integral method Schuh [15] 
calculated unsteady boundary layers on bodies of arbitrary shape and 
for arbitrary variation with time, of the speed outside the boundary 
layer. He observed discontinuity in the slope of the skin friction 
distribution between ~ = 0 and i. 1.0 for the case of a flat plate. 


He attributed this to the failure of the boundary layer theory assumptions 
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near the leading edge. This will have to be resolved by using the 
Navier-Stokes equations in place of the boundary layer equations. With 
the complete Navier-Stokes equations the discontinuity would vanish at 
the leading edge and with it also the discontinuity in the slope of the 
skin friction distribution. 

Schetz and Oh [16] also used the momentum integral technique 
to obtain the velocity field with a new type of profile. They have 
analysed the transient development of the boundary layer on a flat 
plate due to the impulsive start of motion of the surrounding fluid. 

They admit that their approximate solution, valid at every "x" station 
and for all times i.e. 0 < t < ~, tends to mask the detailed mathematical 
behavior of an exact solution. It cannot be expected to resolve the 
problem of "joining" at t = 1. 

Tokuda [17] presented a solution to the Stewartson's problem 
[2,8] without having an essential singularity by "stretching" the co- 
ordinates. His erroneous conclusion that a power series solution exists 
about the Rayleigh solution for small times was corrected by Brown, 
Stewartson and Lam [18]. 

Recently Ban [19] has shown analytically the presence of an 
essential singularity at t = 1. He investigated the velocity and temper- 
ature boundary layers developed on a plane wall by the ideal shock-tube 
flow for the case of weak shock and expansion waves. Following Friedrichs 
theory [20] of symmetric positive linear differential equations Ban has 


provided the proof for the uniqueness and existence of his solutions to 


a } 


edd enizu yd bevipeen on. gt avail Site ataty, built tied. 
AotW . enofteups rave! visbrwod sit to soslq nt enorssups esnos2-rstvel 
ts letnsv’ bivow ytuatsnosetb sft ehottsups eslote-retvell etetqmooserit: 
of? to aqofe ofit nt ituntiaeat arts ozls it dotw bas epbo ontbsel eft: 
notivdiysetb notsatyt able, 
cupindoss fexpstnt muscempmetitubeduoete [Ff] dO bas siete, aatpeten Wt 
sven yenT .sfttovg to sqyd wan 6 Adiw bisit ystoofav ont mbssdo, od) 
tel} 8 no evs! yvebnued sf3 to fnemqofeveb Jnstensyd ond beeyfens 
-btult pnidaucitue sit to notson to Ixste evteluqmt set of oub oe 
nottete "x" yreve te bilsv .nofiufoe atemtxoxggs viens Jeet dtmbs yeaT 
[sottemertem beftetsb edt Azam ot 2bnet .~ > t > 0.9, asmit [1s vot bas, 
— gftt avfoesy of betoeqxs 9d Jonnso +1 .notsuloe 3osxe ns to votvered 
* .f = 7 36 "pntntot" to mefderq 
msfidovq 2'noetiswet2 edt of -nofsuloe 5 betnsesxa [Vf] sbusoT 
-09 aft “patdoserte” yd: ut rs fupnte (sttos2e9 ne patved suodstw [8.S}y 
asetxs notiuloz estvee yewog 5 Jedd notawfonos evosnoyve 2th -2assatbroy 
word yd batoeyies 26w 2omtt [feme rot noftufoe rgtalyeh oft ‘doy 
| [81] me dns nozsrewod2 


— 


= 


née %0 songesiq fit yl feottyisns nwode es [OF] ns8 yldngooh 
~vaqmed: biie-yst sol ovata <bepaptieavitt-eb.f <2 taxpteedupateteiinens 5 
odut-doonde fsebt ent yd Mew ensiq 5 fio beqofeveb 2vayst yvsbauod e1uds: 
aotibaiia patwoliod. .2evew sofansqxe brs Avon Yesw. TO 9269 of'3 107 wolt 
apd ng8 2notteups Isignevatttb sani! evtsteeqoftamege #0 LoSd. ceed | 
© 0f anoisulo2 eid Y sonsteixs bus zeanaupinu end sori Yoong gad bebiverqe — 


the singular parabolic equations. As the basic character of the dif- 
ferential equations is different for t > 1 than for t < 1, he points 
out the necessity of an essential singularity for the problem to be 
properly posed. Murdock [21] developed an integral technique to solve a 
general class of shock-induced boundary layer problems including Ban's 
problem. 

The difficulty at t = 1 was also experienced by Rodkiewicz 
and Reshotko [22] while dealing with the transient weak interactions. 
Their results for temporal shear stress, heat transfer and induced 


pressures cover the range 0 < Tt < I. 


1.3 Choice of Variables 
The literature Survey pertinent to the present work suggests 
that in formulation of the transient boundary layer problem one of the 


independent variables should be of the type: 


The other independent variable has been chosen to be of the form: 


Hey 
n= pes a = | fed . 
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The present problem has been formulated primarily to obtain 
the transient pressure distribution for the case of strong interaction 


(with non-zero pressure gradients). The weak interaction problem (with 
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initial zero pressure gradient) has been dealt with in references [22] 
and [23]. In the present analysis the equations have been developed 
with suitable distortion of the coordinates, which are especially 

suited for the study of flow with pressure gradients. The following 


hypersonic assumption has been made in the theory [22] 


Seg ge Op Heh te teh & (1.1) 
i.e. the free stream adjusts itself instantaneously to the new conditions. 
In a transient supersonic flow, there are three stages in 
the transient pressure distribution before the establishment of an 
Ackeret pressure distribution on the surface of an airfoil corresponding 
to the new free stream conditions. These three stages are: 
(i) the time interval prior to arrival of the fastest signal 
from the leading edge; 
(11) the time interval between arrival of the fastest signal 
and the slowest signal from the leading edge; 
(iii) and the final steady Ackeret pressure distribution after 


the arrival of the slowest signal from the leading edge. 


With the increase in Mach number the time interval between 
arrival of the fastest signal and the slowest signal at a point on 
the body decreases approximately as 1/Me. Therefore, one may say that 
at very high Mach numbers the equivalent of the Ackeret pressure dis- 
tribution is established (from a potential flow point of view) almost 


instantaneous ly. 
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1.4 Approach to the Strong Interaction Problem 

The time-dependent two-dimensional boundary layer equations 
for a compressible fluid have been taken as the governing equations 
for the problem under consideration. These equations have been reduced 
to a coupled set of third order partial differential equations in two 
independent variables n and t mentioned earlier. Next chapter gives a 
detailed account of the transformations employed for such a reduction. 
These coupled equations have been solved for the hypersonic pressure 


gradient 
opted, 
: ve 


which corresponds to the strong interaction case. The numerical method, 
described in Chapter III, is an implicit finite difference technique 
which is inherently stable. This numerical method has bean extensively 
used and produces results of high accuracy. 

The solutions to the momentum and the energy equations give 
at once the shearing stress and heat transfer at the wall for a specified 
wall temperature. The relevant expressions have been obtained in 
section 4.3 of Chapter IV. In order to obtain the viscous interaction 
induced pressure, the boundary layer solutions of Chapter III have been 
utilized to obtain the displacement thickness. Employing Lighthill's 
concept [24] of a one-dimensional piston pushing into air the computed 
displacement thickness has been used to obtain the "piston velocity" 


and subsequently the strong interaction induced pressure over the plate 


8 
| netdov9. notiosys7nl pnowte snd oF flosarqgA af 
ano 3 » seve! yiebrwad leaotansmib-owd Jusbaaqeb~smis onl : 
anotieuns patavevop oft 26 moet need over brult sldtessiquos iw 
: = —_ ; ° 
bayubs Y Nest 1 2noTt7z6 :] g2ai\T Norra yspien Ja tt Le m3) GONG ois oT 


ows of enorssups fsrtnexsttth (efiveg wSb10 bifdt to Js2 beslquea 6 oF 


6 2avfp vedqsdo txeW .variras bsnorsnsm 7 ON& 7 297 CST TS¥ jisbneqobnt 


‘nottoubsy 6 fave vot bsyolons enottamyotensis oft to soveoos bolisted 
r a * a 
U } ora ave sat vot DSviGe nsec §& ' AU | F yer Quod g8¢ sit ~ 
astbsyp — 
a al 
ch 
™ ~ 
3 -borniam Tho! y ont 262 ftosyedal pnowe sd? oF 2baogestiveo dofiw 
oP > A Stn al te » m r «l r 5} | ' 4“. .07 5 “ bea ‘car 
gsuptnfast sang isttib sitnlt tforlqmt ms 2f . Al yesgendo nt badiyoesb 
a ti 
vievfenetxs need 26d bordism lsotvemun 2taT .efdete yltnsiecnt et dotdaw 
hte 
—. : , ae oat ; ‘teas 
Vopywaos tpid to 2tfvee 2osubore- bas beep 
. 4 . = 
evfp enoriéuns yoce adt- bap mutnomom say oF. 2hotguloe snht a 
\ fre) ¥ > ., 


bertissq2 5 10? (low ond 26 Veteneis Jest brs ezeris2 ontveed2 eis sono 36 
é : A AA 


at bantsetdo need svgd enoleesiqxe Jnavelsy snl Mudsnsqnes ‘Tisw 


nord2s1siat eyooehy ant ntatdo oy ees VI a to £.8 noisoez i 
nad sven Li! yotqss2 to 2aofsutoe yavys!l -yisbouod ont essuee9ng deoubat 
wT inset aniyor ong .2eando tnd togmaos! garb ai WIo | best i ttu 


| the 9: bob pridevg petits Jeqotensitib-9ne & to [as] tqsone on 
me S mg ay, : . Agi 2b Hie FLARE, 
veka" ind; Rais at : oh Ts P os a» teh 
a an , ee 
dus ath: [nes S & isis ul Sever peat 


¥ i ee 
, : 


i. 
b 


a 
Ts 


. Pr 


surface. The results for different values of Mach Number, wall heat 
transfer and the transient parameter &(= t/t+1) have been displayed 

graphically for y = 1.4 and y = 1.67. The discussion of the results 
is included in the last chapter. 

It should be mentioned that Lighthill's piston theory is 
accurate only for Mo < 1, where o is the local inclination of the shock 
to the free stream. The changes in the entropy increase with Mo» and 
it becomes necessary to account for them in an adequate formulation for 
large Mo: We neglect the shock curvature and assume the flow to be 
isentropic except across the shock in order to avoid variations of 
entropy across the stream lines. Following Miles [25] we have used the 
extension of Lighthill's piston theory for large M,o- Miles has sug- 
gested that such an extension can be made simply by replacing the Light- 
hill's simple wave relations by the corresponding relations for an oblique 
shock. Expression (4.6) of Chapter IV represents the results of Mile's 
extension as applied to present work. 

Generally the boundary layer in hypersonic flow is thick and 
doubts may arise regarding validity of the Prandtl's theory. Shen's 
[26 ] criterion of 6*/x << 1 has been met in order to justify the use of 


Prandtl's boundary-layer theory. 
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CHAPTER II 
FLOW GOVERNING EQUATIONS 


2.1 Boundary-Layer Equations 

The unsteady, two-dimensional flow of a compressible viscous 
fluid along a semi-infinite flat plate will be considered. With the 
coordinate system fixed with reference to the plate and origin at the 
leading edge, the usual boundary layer approximation gives the following 
equations [29]: 


Conservation of Momentum 
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Conservation of Energy 
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Conservation of Mass 
geri = 2 ard 
22 + 2 (pu) + ar (ov) = 0 (2.4) 
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Qutside the boundary layer where the effects of viscosity and 
heat conduction are neglected, velocity, pressure, temperature and density 
are functions of x and t only. The following relations are satisfied in 


the free-stream: 


OU OU, oU dp 


e . e e 
pipe caer Ce arin ae (2. 1a) 

ah oh dp dp 
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Gy; uj ve aX ) ot Ue aX (2.3a) 


U ) = 0 (2.4a) 
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2.2 Transformed Boundary-Layer Equations 

The equations of Section 2.1 suffer from several disadvantages: 
for example, they may be singular at x = 0, and the boundary layer 
thickness ate greatly with x near the leading edge (in the strong 
interaction regime). In this section they are transformed to a more 
convenient coordinate system by making a change of variables in the 


original boundary layer equations. Some of these new equations exhibit 
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much more clearly the role of certain parameters and of the viscosity 
law, as well as the relation to the incompressible case. They are also 


more suitable for numerical computations. 


2.2a Dorodnitsyn-Howarth Transformation (in the restricted form) 

In the case of zero pressure gradient, Dorodnitsyn-Howarth's 
restricted transformation is equivalent to the Howarth-Stewartson trans- 
formation to be introduced for non-zero pressure gradients. The purpose 
of this transformation is to remove the density p from the formal equations, 


and introduces the following new independent variables 
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The derivatives occurring in (2.1) through (2.4) become 
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The equations (2.8 are now applied to (2.1) and (2.3) to give, 


respectively, 
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one ee ae ee OO (2.9) 
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(2.10) being derived from (2.1) and (2.3) for the new dependent variable 


H with the hypersonic assumption (1.1). 


It is now convenient to introduce the stream function w such that 
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The continuity equation (2.4) is automatically satisfied by this defi- 


nition of the stream function. 


The momentum and the energy equations now become 
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2.2b Howarth-Stewartson Transformation 


For the sake of convenience this transformation will be intro- 


duced in two steps. 


First Step 


The independent variables in this transformation are 


Rete os ty ey (2.15) 
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Second Step 


The independent variables for the second step are: 


Po 


(e) 


Jc Ba Se ete (2.19) 

fe) 
where C is a proportionality factor in the linear viscosity law given by 
C= 0(x) = oT (2.20) 


Expression (2.20) is of the form taken by Chapman and Rubesin 
[30], except that the reference conditions (uo: Te are free stream 
stagnation values. The proportionality factor C serves to match the vis- 
cosity with the more exact Sutherland value at a desired station. If 
this station is taken to be the plate surface, assumed to be kept at 


constant temperature, the result is 


-——) (2-21) 
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In view of (2.19) the following transformation equations may 


be used for the derivatives: 
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with the equation of state for a perfect gas 


Ba ORT (2.25) 


2.2c Final Transformation 


Distorting the time-variable, the independent variables for 


the final transformation are 


ie 

wae 2 ee 

T=) c(2)*(8) df; 275727 (2.26) 
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egal (2.27) 


Use of the above transformation equations for the derivatives in (2.23) 


and (2.24) yields: 
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and the hypersonic assumption (1.1). 


(2.28) 
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(2.30) 


(2,31) 
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(2.38) 


(2.34) 


Equations (2.28) and (2.29) represent the governing equations 


in terms of the new independent variables T, X and Y which are related 


to the physical variables t, x and y by the following expression 
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2.2d Reduction to Similarity Form 
The solution of equations (2.28) and (2.29) could be attempted 
by a number of different approaches. The approach chosen in this analysis 
is to use a similarity transformation. This type of transformation 
changes the two equations into a coupled set of partial differential 
equations in two independent variables n and t. The velocity and the 
total enthalpy profiles at different plate surface locations and different 
times are the same as at any other location and time except for a scale 
factor. 


We assume the following relations for this problem 
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Y = Bx Voan 
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where A, B, C, a, b, Cc, Pp; gq and r are undetermined constants. 


When equations (2.36) are used in the form 
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Ve OVE V 9 ; 
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where the following relations have been employed: 
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Let 
V2 ex" (2.42) 
and B = 2m_ (2.43) 
m+i * : 


Using (2.42) and (2.43) in (2.39) and (2.40) we obtain 
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in (2.44) and (2.45) we obtain after some algebraic manipulations 
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We may finally write (2.47) and (2.48) in the following form 
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For steady state the above equations may be simplified to: 
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Equations (2.51) and (2.52) have been obtained by Cohen and 
Reshotko [31] in a slightly different form. 
Since ue may, in general, depend on X the terms in equations 
(2.49) and (2.50) containing the factor in curly brackets are not yet 
functionally consistent with the rest of the terms in these two equations 
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terms may be achieved by one of the following ways: 
(1) The external Mach number may be a constant greater than zero. 
(2) The external Mach number may be zero (the viscous dissipation 
and compressive-work terms are omitted in equation (2.3) for 
this case). 
(3) The ratio of specific heats y may equal 1 (for most gases, 
this assumption is physically unreasonable). 


(4) The factor 


may be approximately 2 corresponding to hypersonic flow. 


It is the last assumption which would be reasonable for the 
present analysis. The treatment of hypersonic flow also requires the 
introduction of the effects of displacement thickness upon pressure 
gradient. For the flat plate, Lees [32] has shown that the induced 
hypersonic pressure gradient for the strong interaction case corres- 


ponds to 
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As a first order approximation we shall use (2.53) in our subsequent 


analysis. 
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of equations (2.49) and (2.50) alongwith the relation (2.53) we obtain: 
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Letting Pr=1 further simplifies the energy equation (2.55) 
to 
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2.3 Correlation with Reshotko-Rodkiewicz Equations [22] for Zero 
Pressure Gradient 
Equations (2.49) and (2.50), when specialized for zero pressure 


gradient, yield, respectively, 
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For zero pressure gradient "e" condition may be taken as 
reference*, since the local "external" values are constant along the 
*For non-zero pressure gradients, we have used free stream stagnation 


values as the reference conditions, since in the presence of pressure 
gradient the "e" values are not constant. 
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In view of (2.57) the last equation may be further simplified 
to 


2Pr g[l-rf'] - g''-Pr g'[f-2cf] - Pr(y-1) m,¢f'* = 0 (2.61) 


Equations (2.57) and (2.61) have been obtained in reference 
Merale 
From (2.60) we may also write down the following relations 
of interest: 
(a) Static temperature T referred to the free stream stagnation 
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CHAPTER III 
SOLUTION TO THE GOVERNING EQUATIONS 


3.1 General Approach 

The problem formulated in Chapter II requires the solution 
of a third order partial differential equation(2.54) which is coupled 
to a second order partial differential equation (2.56). Equations (2.54) 
and (2.56) are the pertinent equations for the strong interaction case. 
As a closed form solution of these two coupled differential equations 


could not be found, a suitable numerical method had to be adopted. 


3.2 Introduction to the Numerical Method 

The method of solution adopted here, has been developed by 
Clutter and Smith [33] for the solution of compressible laminar boundary 
layer equations with transverse curvature effect. In this method of 
solution the partial derivatives with respect to the modified similarity 
variable €, associated with time and space,are replaced by finite differ- 
ences. The derivatives with respect to the similarity variable n, associ- 
ated with space, are retained. In this way the partial differential 
equations become approximated by ordinary differential equations at any 
particular value of &. 

The particular choice of this method is based on the following 


reasons: the requirement that the method should solve both accurately 
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and rapidly any problem for which the boundary layer equations are valid, 
seems to call for an implicit finite difference technique which is in- 
herently stable and is also exact in the limit. The method chosen here 
comes in this category and has the additional advantage of reducing the 
equation to an ordinary differential equation. Questions regarding the 
existence of the solution, the nature of the solution, and the error 
propagation are much better understood for ordinary differential equations 
than they are for the conventional finite difference methods used for 
the partial differential equations. A final reason for the choice of 
the Clutter-Smith technique is that it is known to produce accurate 
results, has been well explored and does not involve any kind of line- 
arization* used in reference [22]. 

The round-off errors in the computer program can be reduced 
by substitution of f' = f'-1 in the momentum equation (2.54). The 
reason for this is that in equation (2.54) all terms approach zero as 
nN approaches Ne: both S and fi? approach unity and the round-off error 
is primarily introduced when taking their difference. The substitution 


of 


*See appendix A for linearized solutions. 
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in (2.54) gives 


<aCriy) (pe tag (3.2) 
In terms of the Euler's transformation [30] 
a= <r (3.3) 


we may rewrite equation (3.2) as 


= 
(EY (1-6)? SE + CH 1)? - 83 - ae (1-8) 
os 2 2 o> 3= 
b 14 (Qindanpn2et idedl od baptarten) teste. det vatg (3.4) 
n onoe 0& an° an 9n° 


Introduction of the Euler's transformation improves the convergence 


rate of the numerical method. 
A function S may now be introduced in the energy equation 


(2.56) for similar reasons that F was introduced in the momentum equation. 


S is defined as 


Sr@ gf-a15¢ (3.5) 
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Its substitution, alongwith the Euler's transformation, in the energy 


equation (2.56) gives 


-c) S ceqecy (tly — aca ras 
(1-8) Se ECE(E A) + ae(S ne IP SS 


(4) (1-2)? SE + Hs 1)? - G+ 1) - 4e(1-2) 
= 2 = 2 3= 
- (241) REL ae 2 (Peo es (3.7) 
n 


3.3 Finite-Difference Representation of &-derivatives 

This fundamental idea of replacing the &-derivatives by finite 
differences to approximate the partial differential equation was first 
advanced by Hartree and Womersley [34] and subsequently used by Clutter 
and Smith [33]. It may be noted that all of the &-derivatives appearing 
in the momentum and the energy equations are of the first order only. | 

In replacing the €-derivatives at a point by their finite 
difference equivalents, one may use two-point or more accurate three- 
point finite differences. Whereas usually three-point finite differ- 
ences have been used,at the start of a solution only two points are 


available and the two-point form must be used there. 
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Since the momentum equation and the energy equation are para- 
bolic in €, the problem must be solved by marching in the direction of 
E. When the solution has been obtained at all previous stations up to 
and including Emel? the problem is to find the solution at the new 
station Em" The notation system for the finite-difference approxi- 


mation of €-derivatives is shown in the figures la and 1b. 


n 
A A 
i AE 2 7 os 
Em-2 5 Sma] Em ml 2 3 Aenk 


Figure la Figure 1b 


Notation System for Finite-Difference Representation of &-derivatives 


Using a two-point finite difference approximation in the &- 
direction, the momentum equation (3.7) and the energy equation (3.6) 


may be written, respectively, as 


oF a 
on on = = 
ene) Gat Bape ttt CP Gal A (Sy ntl)3 
af _ oF Estee) 
of m,n m-1,n m,n m-l,n, of 
- - +]) (———_——— _)-(-— =) SS | } 
ote en) rt sn z Sm om Sm m1 ane m,n 


ae ‘ ; » f ‘ z Es, 
/ 


~615q 918 nottsupa ygrons oid. ying cicialbilahiaiae Ht 
to nortosyvib od at pntilowsin yd bavio2 sd teum naldarg anit: .3 mf obfod ‘an a 
od qu anita % 2untverq [fs 38 bontetdo nesd 26ft nottutoe eft non 43 ll 
won oft t6 noriuloe sid batt ot af mefdorg ent 73 enibul ont bne q 
~Fxovqgs sonsiattib-ottntt sit vot mod eye nottsion sAT “ne notsete 

.df bos 61 esp Ft oft ot awone et ‘2ovitsviveb~3 to wai : 


df syupra ; Shia _ ST suptt, “int? Be - 
eovisevinab-3 To noftsdasesyqaA soneratt td~sstott Tot maseye “are a 
ea ttn 

-3 oft ni nottemtxorggs soneysttib estat? Snhoqrowst 6 eatau soma 
(3.8) nottsupe ypyene srt bas (1.6) eel. nan 


SY 


oz 3= 
~ f a°F 
“at, th pe | | RM 8g (3.8) 
at et an m,n an? m,n 
(1-& ) «oman, [(1-& ) (41) - 4€ (3t +1)] - as 
1 a ee eT wih i m,n ano m,n 
€ ant 
—-30- 1 3 { msn m-1,4n,4 _ 
an on Eee pny maaan a a a by = 0". (3.9) 


When the three-point finite difference approximation is used, 


the equations (3.7) and (3.6) become, respectively, 
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(Sint) 
Equations (3.8) through (3.11) are ordinary differential equations in n 
with the variable quantities f, f', and S at the m-1] and m-2 stations. 
It may be noted that step size A&, which need not be constant, is not 
a primary parameter; instead, €/A& is. The errors in the two-point and 
three-point formulation, respectively, are 

(nme BF agg Semen) Enrica? oF (3.12) 
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which, for the equal stepsize in €-direction, reduce to 
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*The absolute magnitude of the maximum error in the two-point and three- 
point formulation is given in Appendix E. 
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This means that to have the same accuracy in the solution at 
all stations the stepsize at the second station, i.e. m = 2, must be 
suitably reduced. In the present numerical computations the stepsize 
at the second and third stations was taken half the value used at sub- 


sequent stations. 


3.4 Boundary Conditions Associated with the Equations of Section 3.3 
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When n = Ne? S(n, 96) = 1, and 
S(n,.6) = 0 (3.18) 
h 
When n = 0, S(0,&) = ms and 
. e2 
S(O Fe) Sipe asl (3.19) 
e2 


or for the case of no wall heat transfer 


a = 0 (3.20) 

n=0,€ 

H 
When ve Oe, Sey & a and 

e2 

H 
S(n,0) = p@—- 1 (3.21) 

e2 


We now specify [35] H, such that 2 = 0 ; thus the S-distribution 


1 =0 
at €=0 may be obtained from 
a. 8 +n) = 0 
Semen, GE" 


which upon integration gives 


/ 
*., ee eo pant ae le ee 
TE) as a ov 0 Dey eine 
at a to nea shagiutee | 


bne *. = (3,0)2 = mnie 


| he. eae. * a bore | “8 
eT.e 
(er.£) a i me 
; 


pe [few on to S260 ont Yor 


OS .£ “02 ra 
— sont we ee 


7 - 
a 


| naa 4 
bas . a = (0,o)2 .0= 3,0 2 7 nodW 

ae ag j “4 
rS. | pf. t (0, n)2 OS 
( . ) id Se | af A-« a4 ; > e 2 ace iy 2 | é 


ere 


notiudixtetb-e ont audt ; 0 = ss 2 godt owe -H [8€] yrioege wer 
ft 2 ew 


4] 


$(n 0) = (-S,)———______—__ + § (3.21a) 


This expression has been evaluated using the numerical values 
of Frz9 obtained from (3.17). See appendix B for details of the nu- 


merical integration. 


3.5 Simplifications for an Adiabatic Wall 

The momentum equation is, in general, coupled with the energy 
equation. For the case of an adiabatic wall, however, this coupling 
can be removed. 

The solution to the energy equation (2.55) for the case of 


zero heat transfer at the wall and Pr=1 may be written as 


2, fe eee pa (3.22) 
2 e2 2 e2 ; 


For the problem under consideration the external flow has in- 
stantaneously acquired the new free steam conditions whereas the boundary 
layer itself is non-stationary. 


‘We may rewrite (3.22) in terms of the dimensionless enthalpy as 
S = 1 (3.23) 


which is the desired solution for the case of an insulated wall. Using 
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(3.23) in (2.54) we notice that the momentum equation is uncoupled from 
the energy equation (2.56) for this particular case. A solution to the 
uncoupled equation (2.54) is now easily obtained from the Clutter-Smith 
numerical technique. 

The static temperature Ty at any point in the boundary layer 


after the jump in the free stream velocity may be obtained from 
S(n 50) ee OF 


Rewriting the last expression as 


2 T 2 


u 
h_ 4 
hy +—-=h 


] e2 


we find, after some simplifications, 


, 
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ate eat Oe oftt - £17(n,0)3 (3.24) 


3.6 Procedure for Solving the Momentum and the Energy Equations 
Simultaneous ly 
After the €-derivatives in the momentum and energy equations 
are nec Toced by finite differences, so that the partial differential 
equations are approximated by ordinary differential equations, the problem 
of solution is essentially to find the unknown boundary conditions at 


the wall that satisfy the known outer boundary conditions. The procedure 
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for doing this is given in the following sections. The two equations 
(3.6) and (3.7) are coupled and must be solved simultaneously. 

In order to get a first solution of the momentum equation 
we assume S = 0 at all points in the boundary layer. (As pointed out 
earlier S = 0 or S = ] is the solution of the energy equation for zero 
heat transfer at the wall). This gives us f-distribution as a function 
of n and €. Solution to the momentum equation is now used to solve 
the energy equation with the appropriate boundary conditions for finite 
wall heat transfer. With the new solution to the energy equation we go 
back to solve the momentum equation again and obtain new f-distribution. 
This f-distribution may now be used to obtain improved solution to the 
energy equation. This procedure is continued until convergence of the 
solution to the momentum equation is obtained to a specified accuracy. 
The details of the method are given in the following paragraphs. 

First the momentum equation is solved for 2 < m < m* with 
known values of f from station m-1 and S = 0. The values of f and their 
derivatives from this solution are used to solve the energy equation 
for the same range of m. Thus the S-distribution is obtained, which 
may now be used to replace the earlier S-distribution (S=0). The itera- 
tive procedure is continued until convergence of the solution to the 
momentum equation is obtained. 

In the iterative procedure let Q = i, mip denote the solution 


of the momentum equation with accompanying solution of the energy equation 


*No convergence was obtained beyond m*=8(for y=1.4) and beyond m*=7 (for 
y = 1.67). 
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for stations m=2 to m=m*. 


momentum equation, the procedure is as given below: 


(b) 


Q=0, nlp 


The momentum equation is solved for m=2 to m=m* with S=0 
and the values of f and its derivatives from station m-1. 
At each station m, it is solved by the cut-and-try and 

interpolating procedure to be described in next section. 
The solution may be denoted by fo: 
The solutions fo: the values of S and their derivatives 
from station m-1 are used to solve the energy equation 


for m=2 to m=m*. The solution is denoted by So . 


The method now proceeds to the Q=1 solution. 
2 m* 
Jos i m| 


The momentum equation is solved a second time by using So 
and the solution is denoted by f,. 


Step (b) inQ=0, mp is repeated to obtain S}. 
* 
Oy es lt, m|> 


* 
The procedure in Q = 1, m|5 is repeated, always with the 


latest values of f and S until 


Q= Qhax 
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When this condition is met, the f and S distributions are the 


required solutions of the momentum and the energy equations. Ort was 


chosen to be 6 to give the values of f' from two consecutive iterations 


withine 1077. 


Sale ena eaten Leen 


At the start of the solution assume S = 0 | 


re 


m* Solve Momentum 
= 0, m| Equation for m=2 to m=m*} > 


0 


Solve Energy 
Equation for m=2 
to m=m* (So) 


mes ia a a le a el 


Solve Momentum 
Equation for m=2 to m=m* 


(F,) 


Solve Energy 
Equation for m=2 
to m=m* (S,) 


* 
Q= 1, m5 


Repeat Q-Steps until Q = ney 


Figure 2 Flow Diagram for Solving Boundary-Layer Equations 


3.7 Method of Solution of Momentum Equation at a Particular pestatinnt 

As explained in section 3.3, the replacement of the €-derivatives 
by finite differences results in an ordinary differential equation at 
each E-station. Each of these third order non-linear ordinary differ- 
ential equations are to be solved step by step as the solution proceeds 


in the E-direction. In this method of solution equation (3.7) is solved 
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as an initial value problem with trial values of fu as a third boundary 
condition. A search is then made through the possible values of fF) ‘to 
find the one that satisfies the outer boundary condition that f' approaches 
zero asymptotically as n approaches ot: Before describing the procedure 
for making the search of fi the solution of the equation as an initial 
value problem is considered. It may first be written that 


n 

ria = Oo (git Fis 

Foe {BF ) dn +h (3.25) 
0 


where = (f'') is obtained from equation (3.7) and fis found by the 


searching method. The details of the method of integration to be used 


in (3.25) will be given in section 3.9. The following expressions 


_tria}mso! aden 
f' = | FiMdnes al ict = | f'dn (3.26a) 
re) O 


(3.26b) 


give other quantities required in the solution of (3.7). 
The steps for searching the correct value of on are given 
below: 
(i) A trial value of ae (it was the value at the previous 
station for this particular problem) is used to integrate 


outward from the wall and a check is made to find whether 
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is repeated to yield a high and a low value of ae 

(iii) With these high and low values of es the bounds on the 
correct value of i are further narrowed by splitting the 
difference between the upper and lower bounds and seeking 
the solution again. 

(iv) This splitting-the-difference procedure is used until three 
solutions are obtained such that Fi at n, is between the 
bounds of - eect Sine < k . At least one of the three solu- 


tions must be high and one must be low. 


A three-point interpolation procedure, described below, is 
then used to determine the correct solution that satisfies the outer 
boundary condition F'(n,) a Q, 

If the three trial solutions are 
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1(n) + Apfa(n) + Agf3(n) (3.27) 


where Lagrangian three-point interpolation is used to determine the 


solution that satisfies the outer boundary condition Tn) = 0. Similar 
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relations can be written for f', f'' and f''' with the coefficients Ay» 


A, and A, given by 


A Sie ee eee Spel: Aeere Fee 
{F3(n,) - O(n.) HF (n,) - #3 (n,)) 
cl Fe boo Ge 
A, = Big see i (3.28) 


Accuracy of the solution can be improved by restricting the 
values of the bounds k. 

It may be mentioned in passing that the Lagrangian technique 
is oscillatory in nature and must be used with care: for the problem 
treated here the smooth nature of the function permits its discriminatory 


use. 


3.8 Method of Solution of Energy Equation at a Particular €-Station 
The following procedure is adopted for the inetgration of 
the energy equation. S'' is obtained from (3.6). Then S' can be 


determined from 
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Sy is found by a search technique similar to the one used for solving 
the momentum equation. (S is supposed to be specified here). 


Now S' can be integrated to find S 
S= | S'dn + § (3.30) 


With sy specified, the procedure for searching the correct 
value of S' is as outlined below: 
(i) Equation (3.6) is first solved by using a trial value of cH, 
(which may be taken to be the value at the previous station). 
This solution is stored as $y (n). 
(ii) Depending on whether $y (n,) is greater than zero or less than 


zero, a lower or higher value, respectively, of Sy is tried, 


The second solution is stored as S5(n). 


The two solutions can now be added to give the most general 
solution since the energy equation (3.6) is linear in S. The general 


solution, which can be made to satisfy the outer boundary condition, is 


S$(n) = A S)(n) + B S,(n) 33.) 


en ee Aah 
/ : | P) Rae 

‘antvioe ot bee sno sd of sl tnte ¢ dovse2 6 yd) bawot 
(sve batttosge ad of beeogque et a ne 


| ; tw ees En 


t ya ee im wes | Al 


e] 


| | ie m tg 10 

(0€.€) Frat ) ot oe: Oy), earns 
es = ; ‘ c ; = 

Jasi103 sit pntroisse vot syubss0%q ant ,betttosge Pa ag tw aed: 


‘woled banifuo ee at (2 to suligy 

we to sufev fefyd 8 entey yd bavioe Sart® at (a. -€) noftsupa (fF) i | 
.(nottst2 evotverq aft ts sulsv sit ad ot nadat ‘ed vem lotriw) 
{m5 2 26 beyote at notsutoe efit 

nent eesf yo oves asd? vetbeyp at { jp) 2 verserw ao ontbmeaso (i 


-botve ef ie to Ulovitasqan su ley vedpia to) vomit’ che 
{tr} g2 26 bevote 2i nofiuloe brioose sifT 
: a y aia 
SP A ial 

fsvensg onT .2 nt yeontt at (d.€): or om old gonte p 
et ,mottibnoo yisbnwod ysiuo sit ytertse ot sbem ad mea Ao AW we 
| | rh fu nit " - 
("6.8) E04 (oa mE — 


eo 


a 


farvaitee ~— aid av fp ot babs od won f69., 


Nort fbAo> 


50 


is used along with the condition 
S( Giga sas S,(0) + B S, (0) 


at n = 0 and 


SLO S,(0) = S,(0) 
equation (3.31) gives 
- S,(n,) 
it ees eee (3.32a) 
Sy(n,) - S5(n,) 
with 
Bree eet boeeen) 


In order to obtain more accurate results a bound kK, on S(n,)s 


similar to k on F'(n,)> has been used. 


3.9 Details of the Method of Integration 

The overall procedure for the solution of the two coupled 
equations (3.6) and (3.7) has been described in the previous sections. 
This section is intended to give the details of integration of these 
equations. With the replacement of the €-derivatives by finite differ- 


ences, the problem of solution is essetnially one of integration of a 
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set of ordinary differential equations. The method of integration is a 
predictor-corrector multistep method that uses the Falkner multiple- 
integration extrapolation formulae and the Adams-type multiple-integration 
interpolation formulae, which are described on pages 116-131 of reference 
[36]. The four point form of these formulae is used. The multistep 
method used here requires considerably less computation compared with 
the one-step Runge-Kutta method to produce results of comparable accuracy. 

A special procedure is required to start the integration near 
the wall, the details are given in appendix C. We consider the general 
situation (away from the wall) where the equations have been integrated 
up to n, and we are interested to obtain the values of f and S and their 
derivatives at Mey b= nt An). The extrapolation and interpolation 
formulae are used to approximate the integration indicated in (3.25) 
and (3.26) for the momentum equation. The two-step procedure used in 
the integration is: 

(i) The extrapolation formulae are used first, with the values 
of f''' and f'' at the stations n, n-1, n-2 and n-3, to obtain 
the values of f, f', f'', and f''' at station n+l. The formulae 


_ employed are 
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ntl), = yt Anh, + Age [323% )* 2647! 1 +1597 


1 38F! 3] (3.35) 


where the subscript E denotes "extrapolated". The errors 
in these equations are proportional to (An)? or (An)®. (For 
their exact form, see reference [33])* The Value of fF!" at 
station n+] may now be obtained from the momentum equation 
(3.7) using the extrapolated values of f'', f' and f. This 
is denoted by 


bl 


weauldue )e »F ntl) _? 


! 

n+1). 

(ii) The interpolation formulae may now be used to obtain more 
accurate values of f'', f', f and f''' at the nt] station. 


These formulae are: 


Sua fare oo rid Fit Fius Fius 
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Fags s0F low SPealGFit,, + 19Ft' - SF, + Fi] (3.37) 
n+] eer +1), n n-1 n-2 pg a 
Fi, = + ant + any" [360F!', + 1717!" = 367!) + 7F'',] (3.38) 


n+] 


These interpolated values are finally used in the momentum 


equation (3.7) to obtain the value of f') |, 


urs 4 FCF : Tiel : Fel): 


*A brief listing of these errors is provided in Appendix E. 
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The errors in the interpolation formulae are much smaller in 
their magnitude as compared to those in the extrapolation formulae. 

Also, they are opposite in sign. Therefore the exact value of the vari- 
able in question, say f', must lie within the bounds of the extrapolated 
and interpolated values. This provides a check on the procedure and 

the solution can be made more exact by choosing a small enough stepsize 
rion. 

It may be indicated here that the predictor-corrector technique 
was used only once to obtain the values of f and S at a point in question. 
A more efficient method will be to use predictor once and vary the step- 
size in n such that the corrector is used once or if necessary twice to 
obtain the required accuracy. For the problem under consideration, 
however, it was found that for a five-place accuracy in the values of 
f and S$, the use of corrector was necessary once only. A standard step- 
size of An=0.1] was maintained. 

The formulae for performing the integrations required in the 
energy equation are similar to those for the momentum equation integration 


(refer appendix D for details). 


G2 Olea 7 S| haeeey Mi mt) on re 
/ mie ay 

at vefisme doum ets setumyor nots Toquesnt and 0 sort ae Tha “i 
98 lunrvot nottsToqsytxe ant nf a2orlt of bevsqmoo 26. am Bes | ‘ - 
-tysv st 7 au! 6v do6xe, oft axoTeisdT.. vngke nt st toqgo ons. yadd 02 ; 
betsloqsi3xe eft Yo ebnued sd ntdatw otf geum .'F we .notsesup at ofds | 
bs eywbga0%9 et no YoSr9° 8 2ebtvorq ataT .eaulsy bessfoqvernt i 
esteget2 devons Ifem2 6 pntzoodd yd tasxs ston sbsm sd 162 nottufoe ont 
, ok 
guptniriset sotosvyeo-rosatberq afd Send sven betsotbat ed Yem FT == +k 


* 
i 


om " 


Nortzsup at tntoq & 36 © bas * to zeulsv sad ntstdo od sano yn bea 2 
-qst2 eft wisv bns 9on0 yodotbsyq sau ot ad (ftw borigem tnatottte evom A 
ot sotwi yiseesoen ti 10 gon0 ~b92u et Yososi109 oft dead dowe nm nt este a 
<Moftsyebfenos yebpw me (dovg 9n3 YO .YosIUIDS: bevtupey ons abetdo 

to 2eulsy ont nf Yosyusos soafqeevtt 6 vot Sedd bnuot aew tf eTevewond 
-gat2 bysbnate A .yino sino yrezesoan esw 102307900 Yo oz od? ve brs ¥ r 

| -bantsintem esw f.Oend Yo este 

oi at bevtupsy enottsypatnt ert ontmotveq vot @sTumrro?. onT 
nofdsypesnt no tyeupe musnemom seg 10% geort of el tmta evs nottsupe wrens : 


sem rot aa sone 


ss Ng aan “ert % ‘ r 4 
gow) Me ined mee 


: ’ ¥ 
iy dal 7 
' ’ 


Debris 
oh a hii Wi wy 


54 


CHAPTER IV 
APPLICATION OF THE SOLUTIONS 


4.1 Time-Dependent Displacement Thickness 
The displacement thickness for the time-dependent boundary 


layer may be obtained from the following equation [37] 


> > 
V #6 Ua f (o,U,-eu) dy] 
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co 


+R loa] (eqn) ay] = 0» (4.1) 
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For a two-dimensional boundary layer their result may be written as 
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Expressions (4.3) can be written in terms of the transformed variables 


of Chapter II, namely, 
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Under our hypersonic assumption of aoe 0, the continuity 


equation (2.4a) gives 
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and consequently, expression (4.2) is simplified to 
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4.2 Strong Pressure Interactions 
If we denote the shock angle by o, then for the condition 


that Ma >> al and ous< 1 such .that M,© >> 1, the pressure on the surface 


of the plate may be expressed [25,38] as 
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where Wo is the piston velocity. 


The normal velocity of the displacement surface (piston 
velocity) is given [24,25] by 
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where a,/a,, has been approximated as unity. 
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To be consistent with the steady state results let us define 


an interaction parameter yx by 
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(4.20) 
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Substitution of the pertinent quantities from (4.17) and 


(4.21) in (4.9) finally gives 
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From (4.22) and (4.23) we may obtain the expression for 


the strong interaction induced pressure referenced to the final steady 
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4.3. Shear Stress at the Wall 


The shear stress at the wall may be obtained from 
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Using the’ viscosity law (2.20), we may rewrite (4.30) in the 


following form 
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We now introduce a Reynolds number based on fluid properties 


evaluated at the wall temperature 
iC (4.32) 
and expression (4.31) becomes 
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The advantage of evaluating fluid properties at the wall 
temperature can now be seen from the fact that for the case of 8 = 
the right hand side of equation (4.33) becomes independent of x and — 


X and is a function of-& alone. 


4,4 Heat Transfer at the Wall 


The heat transfer at the wall can be found from 
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or 


oT 
x(—) xq 
rs oy’ WwW 
Nu = pH = (4.36) 
fe) w fe) we oW 


2 Ww e°o W\; Wy @ -] 2 
us SSS SO ie GS) 
Ly ly aM X Yo Po Cp a! i 


rr 
~ eff Us IES (Oo) (ey ICRA 


2 w dinx 
dink -S" (One) 
aTiie “Ty (4537) 


From (4.33) and (4.37) a simple modified Reynold's analogy 


= / Re + ¥ piel 


parameter is obtained 


ee 


R ' 
Cab Sweseat w(ove) ands 
Nu SECCRs ; 
pues 
W 
The numerical results for the transient contributions to 
the shear stress and heat transfer at the wall and to the interaction 


pressure are displayed in appendix F. 
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CHAPTER V 
RESULTS AND CONCLUSIONS 


5.1 Discussion of the Results 

Using the Prandtl boundary-layer equations for a time- 
dependent two-dimensional compressible fluid flow, the strong inter- 
action problem has been analysed. The earlier work of reference [22] 
investigated the weak interaction problem. The numerical method adopted 
there, requires the linearization of the governing equations. For this 
linearization to be valid, the free stream velocity change was restricted 
to about 1%. The numerical method adopted here eliminates this restriction. 
The governing equations, under the suggested transformations of Chapter 
II, are more suitable for analysing the strong interaction case. One 
of the contributions of the present work lies in obtaining equations 
(2.49) and (2.50) which are a coupled set of partial differential 
equations in two independent variables for large Mach numbers. As a 
limit for weak interactions, the present equations reduce to those ob-. 
tained in reference [22]. These coupled equations have been solved for 
B =(y-1)/y, with y = 1.4 and 1.67. The induced pressure results are, 
tyeretore , obtained according to the strong interaction theory. Another 
important contribution of the present work is in obtaining the transient 
strong-interaction induced pressure along with the time-dependent wall 


shear and wall heat transfer for an insulated and non-insulated wall. 
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The distribution of the transient contribution to the velocity 
function fi, is given in Figures 3 through 10 for different values of 
wall heat transfer and for 8 = 0.286 and 0.4. These distributions have 
been obtained corresponding to 1% change in the free stream velocity. 

It may be noticed from these figures that the numerical results have 
been obtained for 0 < € < 0.3 only, since the method of solution becomes 
unstable for — > 0.3. This instability may be explained qualitatively 
if we examine equation (3.7). The leading derivative terms in this 


equation are 
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The above expression has the character of the conduction heat 


equation with the equivalent coefficient of conductivity given by 
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(AyF'+5y+1) and (y+1)/(4yF't5yt1) < & <1, At & = (ytl)/(4yF'+5yt1), K 
becomes infinite. It seems that the discontinuity in the f'-distribution 


occurs upon approaching the neighbourhood of € = (y+1)/(4yf'+5y+1). 
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At the outer edge of the boundary layer, where f' = 0 and the discon- 
tinuity first appears, € has the value of 0.3 for y = 1.4. 

In the discussion of literature in section 1.2 Stewartson's 
Singularity at t=l(or €=0.5) has been described in detail. If the high- 
est order derivative terms in equation (2.49) are considered with B=0, 
the relationship between the singularities at &=0.3 (for 8=0.286) and 
at €=0.5(for8=0) becomes clearer. The highest order and the mixed 


derivative terms are 
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The coefficient of the mixed derivative (5.3) will have different signs 
ror. 0 <1 < Zr and zr < tT < m, Qnce again this coefficient will be 0 
or the equivalent coefficient of conduction will be ~ att =a, . This 
seems to imply that the discontinuity in the solutions will set in for 
t = 1 at the outer edge of the boundary layer where f' = 1. Thus 
t=l1(or —=0.5) is the singularity for 8=0, which corresponds to &§=0.3 
for 8=0.286. | 
The curves of Figure 3 through 10 indicate the beginning of 
a monotonic approach to the final steady state. For 8 = 0.286, Figures 
3 through 6 also show that the approach towards the steady state is 
faster when the wall is kept at a higher temperature as compared to the 
case when the wall is at a lower temperature. (In all these cases the 


flow of heat is always from the fluid to the wall). One possible inter- 


pretation for this trend is related to the effect of wall temperature on 
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the mean density of the fluid within the boundary layer. For the hot 
wall, the boundary-layer density is less than that for the cold wall, 
rendering the hot-wall boundary layer more susceptible to thermodynamic 
changes than the cold-wall boundary layer. This implies that a hot- 
wall boundary layer would approach the steady state faster than the 
cold-wall boundary layer. 

Figures 7 through 10 indicate a similar behavior of (f'-F5) 
for 8 = 0.4. A comparison of these figures with those for 8 = 0.286 
indicates that the approach towards the final steady state is faster 
for 8 = 0.4 for a given value of the wall temperature. 

The transient distribution of the shear function Une is given 
in Figure 11. This figure shows that the approach of the shear function 
towards its steady state value is rapid and monotonic. It may also be 
noted that heating the surface and larger pressure gradients increase 
the wall shear. 

The distribution of the transient contribution to the boundary- 
layer temperature is shown in Figures 12 through 19. The negative values 
in (g,-9) distribution for small n's may be explained if we consider 
an approximate solution to the energy equation (2.56). Such a solu- 


tion for S| = 0 may be written down as 
s(i.t) = S$. + (1-5) f'n.) (5.4a) 


This may be used for Sw ~ Q and small values of 8, since the contri- 


p 


2 


s7 mm rT 
TYO OU) 
’ - » 
- go wt 
1 
3 é 
ie . 
pe 4) 
4 Fel ™ 
Coa. a 
Te 
LE c 
. _% . 
AP ; 
wil 
f 
‘ ‘ 
‘ a? 4 
cw 
> a ie 
i wod 3} wl 
ay r4 of 
SV TIS6PSsii 


reno 


_ S232 5 Ff 
~ 

- 5 

sh .2) 

wi: 


at 
i 


oy) 


7 


oe 


e 
4 


(r,t) '? 


rhe 


is ort? tend ewode siugtt arnt tf rupit at 


3 ‘ . - | = jae 
vows Sf esyuptd nr nwode at yseraged 195 i 


nbiveups. ypiens snd - nor sutee ital 1d G6 16 


hs 02 ans santa .é9 to © eaviay (fame: ‘nb ¢ 
uy “a : | a 


. nan Foal - 
od odd ntdtrwebtult and t0 waren Sem asl 
ned jee ad, ybe argh : ysl pias ‘oft tts Tew 


nit oft 2bvewo! toeoves oid teas 2edsoro 


oyom Yayel ywaedtved -Tfaw<tod ws ontab bas 


= } J : ' ; mI “ a 
bv a! ‘Yrsehhoe Ph ew-p Tes ait 6 id wd et 
E , } ( ¢ it t @ 
ta si3 sasonqgs bison veyel yas ees sw 
“h 


f1QNe PMY sboudd T Isw-bfieo, 


stsothnt Oh apueris: § seawert « — 


a). 
a4 , — ; a et 
S°Hy p rT 925i te wmelit Cid A 6.0 ° a a? 4 


a 


fisw oft to subey ‘aiied 6 wT §.0 = 9 YOT- 
; ° ) , - 
ioOrsvaoTise ‘b share nav aT 


1 biqey 2t ouisy ajed2 ybaste zit 2b16 wold ; 
-f wh v :] 
rgovel bie Sastre. sit entised et it 


t 3 
4-4f r 


Naga Fis er | 


Di uGra garb aAT i iw 


* 3d. Yam Pb [ Laine 4OT nor udtraat | dn 9 ») a 


Wie A dy 


u 


69 
bution of g(f!2-s) to the momentum equation is negligible as compared 


to the other terms. For S,, = 0 we obtain from (5.4a) 
Shs fe (5.4b) 
For large Mach numbers expression (2.60) also gives 


(yng) = SEL? (s,-#52-se 012) (5.5a) 


Using (5.4b) in (5.5a) we obtain, for Sy 90 


(go-g) = OCD) y? (tient odtibt ho lum | (5.5b) 
For the problem under consideration (f'-f5) is always positive. How- 
ever, the expression in the curly bracket is negative for small n's and 
positive for large n's. This gives the qualitative explanation for 

the "peculiar" behavior of (95-9) at small n's. In order to explain 

this physically we consider the viscous dissipation close to the plate. 
The slope of the velocity profile which is a measure of the shear stress, 
changes very fast for small values of n's and &'s. This means that there 
is substantial viscous dissipation for this range of n and &. This raises 
the temperature of the fluid next to the plate higher than the final 
steady state temperature. For large values of n's however, the velocity 
profile is less steep and accordingly the viscous dissipation is not 


very significant.. Therefore, the temperature of that part of the boundary- 
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70 
layer fluid is smaller than the final steady state temperature. As time 
increases (i.e for large values of &'s) the heat transfer in the 
boundary layer helps to bring the temperature of the fluid towards 
the final steady state temperature for all values of n's. Figure 12 
through 19 also indicate that for a given Mach number and wall temper- 
ature, the transient contribution to the temperature in the boundary 
layer is larger for 8 = 0.4 as compared with 8 = 0.286. 

The linearized solutions for the transient contribution to 
the wall shear, boundary-layer temperature and the velocity function 
for the case of an adiabatic wall and y = 1.4 have been presented in 
appendix A. It has been found that the wall shear predicted by the 
linearized solutions is about 1% lower than the one obtained without 
linearization using the Clutter-Smith technique. The values for the 
transient boundary-layer temperature are smaller by about 1/2% using 
the linearized solutions. The agreement between the linearized and 
the non-linearized solutions appears to be very good for the 1% change 
in the plate velocity. For larger changes in the plate velocity the 
linearized solutions are expected to be in more error as compared to- 
the non-linearized solutions. 

The transient distribution of Reynolds analogy parameter 
is given in Figure 20. These curves resemble the f''(0,€) curves 
(Fig. 11) because of the relatively small variation in magnitude of 
S'(0,8)/(1-S,,) compared with that of f''(0,€). Utilizing expression 
(4.38) and Figure 11 we may easily obtain the heat transfer at the 


wall from Figure 20 for any particular case in question. 
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The strong-interaction transient contribution to the pressure 
distribution for the different values of wall temperature and Mach 
numbers is given in Figures 21 through 26. We notice from these figures 
that for a given Mach number and the pressure gradient parameter 8, 
the transient contribution is larger for a hot wall as compared to that 
for a cold wall. We may thus conclude that cooling the surface tends 
to induce a smaller pressure rise in the interaction zone. 

With the increase in Mach number the transient contribution 
goes down for fixed values of S., and 8. For a given Sy and Mach number, 
the transient contribution is smaller for 8 = 0.286 as compared with 
8 = 0.4. 

Some of the induced pressure results have been further 
analysed in Figure 27 for Ma = 5 and xy = 6. This figure indicates 
that the transient contribution decreases faster for large values of 
8 and wall temperature. 

The numerical results pertinent to expression (4.22) reveal 
that for most of the cases under consideration, the following approxi- 


mate expression may be used: 


2 
P=Da (3y-1) M 
ogee ta cal eramegre 
ik Mo eo G6 
ae 
-1)"° M 
en e px]? Ga 


oy! ‘ 7a é 


. ‘ s : a2 ng or ~ - wy . iy 
oid of nottudiyines jrigtens1d notsseyssnt-gnow es aa, 


ul . ‘ 


foam bis ayussraqiss ‘isw to esutsay Inetstt Tb saad 0t, oot nti ‘date 


>a vit 
> Tie - 
- 4 
+ 
ad 
ale. , 
~~ ior ive 
‘ ‘7 i =) 
4 
hia Mi 
© 
2 ¢ : 
€ 


‘ > Ne ! a 
iioyv} SoTsin OW -.4 uoVdh FS esvuery vt riav tp. at eve ei ee 

2 i : - bral > 
naveq Fosibs yD 3 319 3ftt bap jodnud foe e0 rp 6 seh ey S 
‘7 


noo 26 {few tod 5 Yot Werte! 2h ‘SYSENINOD snshanent of - 
VSy) OW ir bfoa om OT 
-9nes floflosieint oAs af. sary Sivezarq val Teme 6 soubs 


nsfens yt eft ssdmun dosh ar sesoyont ans AvTW 1 oe 


bins 2 navig 6 107 «8 DnB. 2 7o aouTey oxtt yo? nwob 4 


mio 26 o3S.0 = 8 X97 ifeme 2h nor judtad nos. toot enens 


7 
cody _ - 
0= 


‘ut Aged sven 2 [pesn SW 25g bsaubaf sft to smoe . Av 


Dnt Syuptt 2rd. .9 = x bas @ = aft YOR - suptt af bevelang 
"éi Yot 193 2st aozct ro9b not dud ryan SnsT air ont 38 ft: 
SE togiEs | Tow t sa 

1X9 GF Ineniivsq etiuess (69 teva wn ent Par 


\ of oF six ii Sseyvysptenos  Sbnu 259 ont. 14. $200 "0 + tend Ag 


i ta 


~sbeeuad ism not eae > 


3 M (f=e) .. ¢ m F B Z 


Peete + of.) Oe oY os 
aM 5 j « gt 


Ta 


5.2 Concluding Remarks 

The present results may be utilized to obtain the history of 
the temperature distribution at the wall, the transient drag force (ob- 
tained through the integration of the shear stresses at various times), 
and the transient induced pressure over the wall. 

In order to extend the present solutions beyond & = (yt+1)/ 
(4yf'+5y+1), a perce approach should be adopted. Ban's [19] method 
of solution appears to be more plausible for this region. Equation 
(2.54) although parabolic in nature, may be subjected to the "elliptic" 
boundary conditions in the region (y+1)/(4yf't+5y+1) < & <1. The 
present solutions evaluated at € = 0.3 (for y = 1.4) may be used to 
provide one set of boundary conditions. 

Once the solutions for the above range of € have been com- 
pleted it would be desirable to extent the analysis to cases of Prandtl 
number other than unity. Similar investigations for other airfoil 
profiles like an inclined flat plate, circular cone etc. should also 


be of interest. 
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APPENDIX A 
LINEARIZED SOLUTIONS FOR AN INSULATED SURFACE 


A.1 Linearized Form of the Problem 


Let 
F(n,t) = f,(n) + Af(n,t) (A.1) 
where fo satisfies the steady state momentum equation, namely, 


fi + FyfLY = B(FS+S5) = 0 (A.2) 


2 Ge 


The momentum equation (2.54), after omitting small terms of the second 


order, becomes 


sat tatty") 


(CY att + 204) tyart acl : 


Y 2 
Vt 4 wt Lael y-1 12 a blag ta eet bglels yl & 
- fof = fylafie apts + (Leyes = fyfyt = #ytt= (Ys = 0 (A.3) 


However, if we prescribe H, so that S(n,0) = 1 then [60] 
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satisfies the energy equation (2.56) with its space conditions, and 


consequently the linearized momentum equation becomes 


vaas | hel ! i res. ee a vt 
( * AT eet i ) foAf At (fAAf fh Af) foAf 


ar nee ete (A.5) 


with its steady state counterpart 
gia yas nt ach ea (A.6) 
2 22 Y Z : 


In terms of the Euler's transformation (3.3) we may rewrite 


expression (A.5) as: 


+1) a(Af') “ly ging » O(AF') _ 6 BC AF) 
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Paton fo At -AATe nD (A.7) 


The boundary conditions associated with equation (A.7) are 


(AF) (ng.é) = 0 (A.8) 
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(Af) (0,€) = 0 (A.10) 
(AT) “tne .0} = €: “ed (ar. “e0) = f! (ACT) 
0 Ve2 Ue2 ¢ 


A.2 Numerical Solution to the Linearized Equation 

The solution to equation (A.7) with the boundary conditions 
(A.8) through (A.11) has been obtained using the numerical method de- 
scribed by Rodkiewicz and Reshotko [22] with certain modifications in 
the computational molecule. In the numerical computations of reference 
[22] the 7-point computational molecule, indicated in Figure A.1, was 
used for all values of E(ort). In order to obtain more accurate results 
in the present work the 7-point computational molecule was used only 
for the first step in &-direction. For the subsequent values the 10- 
point computational molecule indicated in Figure A.2 was used with 


variable stepsize in &-direction. 


3 


Figure A.1 The 7-point Computational Figure A.2 The 10-point Computational 
Molecule; H = 0.05, Ty = 0.025 Molecule; H = 0.05, T = 0.05, 
Ty = 0.025 
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The momentum equation (A.7) may be written in the finite 


difference approximation as: 


ann Af ntl : Den At nan cman Af n= dian At n-2 r emyn (A. 12) 
where for the case of 7-point computational molecule we have: 
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While for the case of the 10-point computational molecule we get 
Sy eed eee : i 2 ' 
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In the foregoing the following finite difference approximations have 


been used, 


For 7-Point Computational Molecule: 
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: 3 
Expressions for = . : af and eC! are similar to those for the 
on on 


7-point computational molecule. 


Equation (A.12) may also be written in the form 


Af = Rae Af 


m,n m,n mnt] i. Ch (A.30) 


aa 


which was obtained by considering the boundary conditions (A.9) and 


(A.10) with Rath and Ch. defined respectively as: 


=F a 
msn 


Ran . (b an ale d R ) (A.31) 


+ 
m,n m,n m,n-1 myn m.n-1'm.n=2 


- d 
C 2 enn ene 
+ 
mn | Den CH, 


n®mn-2¢m,n-1" 
n¥nn-18m sn-2 


m m, 
fowshel om 


where 


Letting n = N in the free stream (where the 'e' condition 


prevails), we may write, from the boundary condition (A.8), the following 
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Ce Tim (> 3AF, ttf n1-4fa nz) = 9 (A.33) 
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Further using expression (A.30) for Af a N-1 29d AF, yo in 


relation (A.33), we obtain 
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Cn N-2 


C a 
m,N-1 4-R 
Af aN = IE (A.34) 


4-Re N-2 m,N-1 


Expression (A.34) may be utilized at the boundary-layer edge (i.e. at 
n=N) as the starting point in the use of recurrence relation (A.30). 

In actual computations the following procedure is used. First 
we calculate the coefficients Rien and Gren? It may be noted that ex- 
pressions (A.31) and (A.32), when specialized to station m, require 
the knowledge of the values of function Af at station m-1. The Af- 
distribution for zero-time(m=1) is known from expression (A.11) and 
consequently we may calculate the coefficients Ran and Cnn for the 
Station m=2. Once these are found for all values of n we may obtain 


Af , from expression (A.30) by incorporating expression (A.34) at the 


2, 
outer edge. This procedure was repeated at other values of m covering 
the range 0 < —€ < 0.3 (for y = 1.4). For & > 0.3, the method became 
unstable. 

The Af-solutions were used to obtain the transient contribu- 
tion to the velocity function (Af'). The distribution of (Af') is 
shown in Figure A.3. It may be noted that each curve has an inflection 
point which penetrates into the region of higher values of n as &€ in- 
creases. 

The time-dependent shearing stress at the wall is presented 


in Figure A.4. Figure A.5 gives the transient contribution to the 


boundary-layer temperature for various values of & and for i = 5, 8 
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Figure A.4 Distribution of the transient contribution to 
the shear-stress for 8 = 0.286 ana no wall heat transfer 
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Figure A.5 Distribution of the transient contribution to the 


temperature distribution for 8 = 0.286 and no wall heat transfer 
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and 10. 

The linearization used to solve the momentum equation re- 
quires that the change in the free stream velocity should not be too 
large (say, about 1 to 1.5%). A comparison between the linearized 


and non-linearized solutions is indicated in the last chapter. 
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APPENDIX B 
NUMERICAL INTEGRATION TO SPECIFY INITIAL 
CONDITION ON TOTAL ENTHALPY 


B.1 Details of the Numerical Integration 
In order to evaluate the expression 
- ae: )d 
=) (f-2ntny an 
| eo env dn 
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e I (Fragtn)dn 
fe 
fe) 


dn 


with the numerical values of obtained from (3.17), we write 


E=0 

f Feeoin)d 
ey a + E=0 n)an 
$'(n,0) = $'(0,0) e 

= $'(0,0) Feag(n) 
ee 
where 3 (F,.gtn)dn 
Frain) =e. 


From (B.2) we may also write 
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00 0 ee eee (B.4) 


Now, 


may be evaluated numerically by the trapezoidal integration formula. 


Therefore, we write 
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Here An denotes the stepsize (=n(n) - n(n-1)) in the n-direction. 
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Subscript 


E=0 has been dropped from the function F(n) for the sake of convenience. 


Furthermore, we may write 


n 
Sn.) = { 5'(n,0) an + 5, 
0 
or Nn“ An 
S(N*,0) = 5,_p(N*) = | Sig ants, 


0) 


Using the trapezoidal rule for integration we obtain 
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Expressions (B.9) may be used to evaluate expression (B.8b) which gives 
S(n,0) distribution represented by expression (B.1). This is the initial 


condition used on S$ in the numerical solution. 
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APPENDIX C 


STARTING PROCEDURE FOR THE INTEGRATION 
OF THE GOVERNING DIFFERENTIAL EQUATIONS 


C.1 Starting the Solution near the Wall 


The extrapolation-interpolation formulae in section 3.9 of 


93 


Chapter III require values of the variables at four previous n-stations. 


To start the integration at the wall Taylor's series has been employed. 


For f'' we may write 


Fin ie ainda rains (an) iv eye (C.1) 
In the present computations, only two terms of this expansion 
have been used. In order to have the same accuracy as in the extrapo- 
lation and interpolation formulae, however, much shorter steps in n are 
to be used in the Taylor's series. The study made in reference [33] 
shows that a five-place accuracy in the values of f and S can be main-- 
tained by using the Taylor's series to obtain these values at only 
n= a from the wall. Here An denotes the stepsize used in the extra- 
polation-interpolation formulae. The numerical value of An used in 
the computations is 0.1. The two-point and three-point extrapolation 
formulae are used to build the values up to the full length step size 


An. Again the accuracy requirements call for step sizes of = and se 
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to be used for the two-point and three-point forms respectively. 


The 


expressions to be used with different step sizes are summarized below. 


(i) Immediately next to the wall (i.e. at n = An/16) Taylor's 


series is used with the step size An/16: 
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(ii) At n = An/8, 2-point extrapolation formulae are used with 


the step size An/16: 
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(iii) At n = An/4, 2-point extrapolation formulae are employed 
with the step size An/8: 
Fans 7 Tange * Te Faye tw He 
Fan/A = fing t Te [3F yt yee fh] (G.. 
Fan/4 = Tanya t B Tine * + Loy [af ne tw (C. 
Sina = Sinsa * 18 ESinye7Sw') (C 
Sansa = Sania * 18 Sansa Sw (C. 
(iv) At n = An/2, n -3 An and n = An, 2-point extrapolation 
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(v) At ne 3 an, n= 2An, n= 2 An and n = 3An, 3-point 


extrapolation formulae are used with the step size An/2: 
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For n > 4An, the 4-point extrapolation-interpolation formulae of 
section (3.9) are used with the regular step size of An. 

Use of Taylor's series to start the integration near the wall 
has the advantage that it can be checked in a simple and effective 


manner for accuracy. If Runge-Kutta integration technique were used 
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as a starting method, the error estimation and checking for accuracy 


would be difficult. 
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APPENDIX D 
INTEGRATION FORMULAE FOR THE ENERGY EQUATION 


D.1 Extrapolation-Interpolation Formulae 

The 4-point extrapolation-interpolation formulae for carrying 
out the integrations required in solution of the energy equation are 
similar to those required for momentum equation, described in section 


3.9 of Chapter III. These formulae are: 


Stl). = 51 + SF [55 51'-505'' +375! '-95'' 3] (D.1) 
Sn), = 5 + SF [555'-595'_1+375"_,-95"_4] (D.2) 
and from the energy equation (3.6) 
ee irs F(S iat) » S41) ) 
E E E 


where the subscript E denotes "extrapolated". 


The interpolated values are obtained from 


Gi oo el An vt Cll felt set.) : 
Snt1 = Sn * 34 9541) 98 55217 Sn-2! (D.3) 
55 + MDs | 
S41 = Sy, +35 [9S +198) -55) 4* 1 Si 2d (D.4) 
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APPENDIX E 
ERROR INVOLVED IN NUMERICAL INTEGRATION 


E.1 Error in Finite-Difference Representation of &-derivatives 

The errors in the two-point and three-point finite difference 
approximation of €-derivatives have been given in Section 3.3. These 
are indicated by expressions (3.12) and (3.13). The maximum error 


involved from this source is given below: 


2 
Two Point Formulation: aes = 0.00085 (E.1) 
o& max 
(AB\e gaek 
Three Point Formulation: ae = 0.00075 (E.2) 
3 9e9 max 


E.2 Error in Integration Expressions 
E.2.1 Extrapolation Formulae 
The errors involved in expressions for f'', f', and f, re- 


spectively, are: 
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E.2.2 Interpolation Formulae 
The errors involved are 


In expression for f'!': 
19 Fu 
Eg < - 35g (an) °C '(n) J 


In expression for f': 


In expression for f: 


Eg < - gage (an)? Fn) 


E.3 Accuracy of the Numerical Results 


(E.6) 


(Ee dw) 


(E.8) 


The quantities that affect the accuracy of the numerical 


results are An, Ne? k, A& and tae Since Cohen and Reshotko's [31] 


results for the steady state similar flows are the only highly accu- 


rate solutions known, these were first used to check the programming 


of the equations. The four-place accuracy of Cohen and Reshotko's 
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solutions was obtained by choosing the following values for the above 


mentioned parameters: 


Ages O02 lens = 5, ke= 120 (2/ Ae) i205 


e 


and 


fol 


Ty i 1 
: + a ot “* £. 
' 1% st - 
— t 
Patt 3 4 Ae 
“908 VIAPIA VINO Shs Si 
) | j 
S07q ON. Jos 
= o te - 
: 129% im | 
~~ + s n ‘ 
I 7 
ig 4 ' | < wi ’ 


yee y> (aa) 


ct oeee: tf ; ye 
ag A Mt ee See 7 D t 
so" “(pid) oaT ~ 2 33 


ho Yoeawons: 


Vi . mabe: 
i Lae é bas * 
: r 9 : : 
7 hy ! 


ahi Fiero? sotsefora 


: 
1 ge se aa 


a1 bavfovit evor"s 


t 
{ oun 
ap co my ay —— 
1 "F yo? noteeetqxe al 


a . 
+f my ; 
TY i ~~ at 
\ “| pA Pp s 
j 7 aA) Sy -) a4 ; 


| sani Gh 
sotramiAt-oft to yoswook £.9 

; Oe 
Foote. tend esists snBup sat ve 


a” oft on aif ves 


2ifuweon Ff 


Pa bas 3A ets. 
2 atete, ybsote ort not ash ae 


ba 2 iy J2 “f 7 si ‘9M g2ond cava senot uta 2 90s 


us tiel 


ae Yuet ail Sonal ber oe | 


== 


102 


APPENDIX F 
CURVES FOR CHAPTERS III AND IV 
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Transient Distribution of Shear Function 
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Figure 12 Distribution of the Transient Contribution to the 
Boundary-Layer Temperature for Sy = 0 and B = 0.286 
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Figure 13 Distribution of the Transient Contribution to the 
Boundary-Layer Temperature for Sy = 0.2 and B = 0.286 
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Figure 14 Distribution of the Transient Contribution to the 
Boundary-Layer Temperature for Sw = 0.6 and B = 0.286 
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Figure 15 Distribution of the Teeter t Contribution to the 
Boundary-Layer Temperature for an Insulated Plate and 8 = 0.286 
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Figure 16 Distribution of the Thransient Contribution to the 
Boundary-Layer Temperature for Sy = 0 and B = 0.4 
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Figure 17 Distribution of the Transient Contribution to the 


Boundary-Layer Temperature for Sy =O 2eand 6a 014 
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Figure 18 - Distribution of the Transient Contribution to the 
Boundary-Layer Temperature for Sw = 0.6 and B = 0.4 
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Figure 19 Distribution of the Transient Contribution to the 
Boundary-Layer Temperature for an Insulated Plate and g = 0.4 
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Figure 20 Transient Distribution of Reynolds Analogy Parameter 
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Figure 21 Transient Contribution to the Induced Pressure 
Distribution for Sy = 0 and M, = 5 
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Figure 22 Transient Contribution to the Induced Pressure 
Distribution for Sie! and M, = 10 
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Figure 23 Transient Contribution to the Induced Pressure 
Distribution for S., = 0.6 and M, = 5 
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Figure 24 Transient Contribution to the Induced Pressure 
Distribution for S = 0.6 and M, = 10 
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Figure 25 Transient Contribution to the Induced Pressure 
Distribution for an Insulated Plate and M, = 5 
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Figure 26 Transient Contribution to the Induced Pressure 
Distribution for an Insulated Plate and Ma = 10 
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Figure 27 Transient Contribution to the Induced 


Pressure Distribution for y = 6 and = 5 
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